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Abstract 



We determine approximate formula for the ground state energy of anyons 
in 2D parabolic well which is valid for the arbitrary anyonic factor v and 
number of particles N in the system. We assume that centre of mass motion 
energy is not excluded from the energy of the system. Formula for ground 
state energy calculated by variational principle contains logarithmic diver- 
gence at small distances between two anyons which is regularized by cut-off 
parameter. By equating this variational formula to the analogous formula of 
Wu near bosonic limit (y ~ 0)we determine the value of the cut-off and thus 
derive the approximate formula for the ground state energy for the any v and 
N. We checked this formula at v = 1, when anyons become fermions, for 
the systems containing two to thirty particles. We find that our approximate 
formula has an accuracy within 6%. It turns out, at the big number N limit 
the ground state energy has square root dependence on factor v. 



I. INTRODUCTION 

Anyons - two-dimensional particles obeying fractional statistics are one of possible quasi- 
particle excitations in fractional quantum Hall effect and High T c superconductors 

The determination of the energy spectrum of anyons in 2D parabolic well is actual and 
interesting task of quantum mechanics. The analytical calculation of the energy spectrum of 
the two anyons in 2D parabolic well was performed in |6|.|7| . The energy spectrum of systems 
containing more than two anyons in the same external field can not be obtained analytically. 
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Wu H generalized the class of exact two anyon solution into three anyons containing system 
and found it's energy spectrum. There are the numerical calculations of the lowest part of 
the energy spectrum of three PJT0|] and four [11] anyons in 2D parabolic well. 



There is separate quantization of centre mass motion energy and relative motion energy 
to each other of anyons in 2D parabolic well The ground state of relative motion 

energy of the systems of two, three and four anyons can be linearly interpolated between 
bosonic and fermionic spectra [p|-|TT|| when anyonic factor v changes from to 1. 

Linear dependence of full ground state energy ( without separation of cenre mass motion 
energy ) from anyonic factor v exists for arbitrary number of anyons only near the bosonic 
limit of spectra |J, i.e. at v — > 0. Beginning with the case of three anyons it deviates from 
linear character near the fermionic limit of spectra ||, i.e. at v — > 1. 

By variational calculation we determine approximate formula for the ground state energy 
of anyons in 2D parabolic well which is valid for the arbitrary anyonic factor v and number 
of particles N in system. We assume that centre of mass motion energy is not excluded 
from energy of system and anyons represented by bosonic variational wave function. As in 
Laughlin's variational calculation of infinite system of free anyons |I2]| , our formula for the 
ground state energy calculated by variational principle contains logarithmic divergence at 
small distances between two anyons regularized by cut-off. We equate our variational ground 
state energy formula to the same formula of Wu || for anyons in 2D parabolic well near 
the bosonic limit on the factor v to determine the cut-off and thus derive the approximate 
formula for the ground state energy at the any factor v and number N. We check this 
formula at v = 1, when anyons become fermions, for the systems containing two to thirty 
particles. We find that our approximate formula has an accuracy within 6%. It turns out, 
at the big number N limit the ground state energy has square root dependence on factor v. 

The following section 2 describes the hamiltonian of the system of anyons in 2D parabolic 
well and the choice of trial wave function for the variational calculation and physical argu- 
ment for this choice. Section 3 contains results of variational calculation. 



II. HAMILTONIAN AND WAVE FUNCTION OF SYSTEM 

The Hamiltonian of N anyons system in 2D parabolic well has a form: 

1 N Mi i 2 N 

zm i=i z i=i 

Here M - mass of particle, p = — iaV, where V = i— — h j— , uj - eigen frequency of 

ox ay 

oscillatory motion of free particles in 2D parabolic well and r*j - radius vector of i - th 
particle number. 

Vector potential for anyons A{fi) m (S) ^ s 



N — -i 

z x r. 



^) = ^EVrf- ( 2 ) 



j>% i r y 



Here z is unit vector perpendicular to 2D plane where situated system,^ = r. \ — fj and 
anyonic factor v characterizes the form of fractional statistics ( spin of the anyon ): it can 
be changed from v — - noninteracting bosons to v — 1 - free fermions. 
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In the variational consideration the ground state energy of system determined as mini- 
mum of energy 



E = J^ T (R)H^ T (R)dR 

Here H - the Hamiltonian, ^/t{R) - trial wave function of system that has a variational 
parameter and R - the coordinates of all particles. Minimum energy of E (|3]) is reached by 
variation of variational parameter of trial wave function. 

We suppose that trial wave function is provided by mean field approximation of Fet- 



ter, Hanna and Laughlin [|13j. In this approximation particles move in the mean field A: 



A(r) = -B x f (4) 
2 

induced by the average density p of particles in the system where B = 2itphvz, p = — - and 

r a - mean distance between particles. According magnetic field B one can introduce cor- 
responding magnetic length an = (h/B) 1 ^ 2 and cyclotron frequency uj c = B/M. Energetic 

spectrum of Hamiltonian (P with A and 2D parabolic well has the same energetic spectrum 
structure as a free particle in 2D parabolic well [14j . 



Thus assuming a bosonic representation of anyons 

N 

* T (R) = II *r(ri) (5) 
i=i 

for the ground state of anyons in 2D parabolic well one can bring wave function \l/r(^) in 
the following form: 

*t(*) = Cexp (-{a + ^) (X ^ 2) ) • (6) 

Here C - normalization constant, a variational parameter and there assumed that mean 
distance between particles r a equals of R a = {h/ Moj ) 1 ^ 2 - raduis localization for the wave 
function in the ground state of 2D oscillator. This equality follows from fact that interaction 
between two anyons at any not zero anyonic factor v has repulsive character and relation 
S n = miRl where S n - area that occupies by free particle at n - th energetic level of 2D 
oscillator. Last relation for the area S n one can easily get by using a virial theorem for 2D 
oscillator. 

It is convenient to express all energetic quantities by energetic quanta of 2D oscillator 
Tiujo and length quantities by r Q . 

Normalized trial wave function of system has a form: 

^t(R) = [—^J nexpf-(a + */) v 2 1- (7) 
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III. VARIATIONAL CONSIDERATION. RESULT OF CALCULATION. 



At the calculation of variational energy E (|3|) it is convenient to deal with quantity 

e l (r) = VT\ii)m T (R). 

As trial wave function ^t{R) is not true of the Hamiltonian ([I]), so 

E L (R) = ReE L (R) + HmE L {R). (8) 
By the acting of the Hamiltonian H ([[]) on ^t(R) (|7D we find 

-/V 2 2/ \ 2 2 

ReE L (R) = J> + v + - + y}) + V - {A{n))\ (9) 

8=1 

and 

v 

Im^i?) = -v{a + v) ^2(A(fi)fi). (10) 

i=l 

In the expression ( |10|) we have a scalar product of vector A(f^) and fj. 
A variational energy ([$]) expressed by El(R) is 

E = J Vt(R) E L (R)V T (R)dR. (11) 

Direct calculation gives 

/ * T (-R) ImE L (R)ty T (R)dR = 0. (12) 

So, on energy E contributes only ReE L (R). 

Integrals ( JTT| ) with ^ T (i2) (0) of first three terms in square brackets (Q) are calculating 
elementary. Problems appear when we are calculating integrals with term proportional of 

(A(fi)) 2 . This term gives two kind integrals: integrals like / ^t(R) ,-, , 9 ^T(R)dR where 



\Ti 3 



1 1 ->. - 

i 7^ j number of which are iV(iV — 1) and integrals like / ^t(R) t-^— ["p - :^T(R)dR where 

1 | | ^ik | 

' 7^ Ji ! V ^ an d 3 7^ k number of which are N(N — 1)(N — 2). 

By taking account of integrals ( see Gradshtein and Ryzik [|15] ): 



Ei{ax)e-^dx = — In ( £ - 1 ) , (13) 



'0 /i \a 

where constants a > 0, Ren > and // > a and 

f° E^-^e'^dx = - In f £ + 1 V (14) 

where constants Re{(3 + /x) > 0,/i > and E^y) = —vp J^e~ lz dz/z - the exponential 
integral where 7 > and vp J indicates integral in main quantity mean one can get 
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J ^ T (R)^y T (R)dR = (a + u) In (?) , (15) 



where cut-off parameter 5 regularize the integral (15|). 
And an other hand the calculation gives 



f V T (R) V T {R)dR = (a + u)A, (16) 



where numerical value A ps* 1. 

So, term by term averaged quantity ReE L (R) (|) gives expression for energy 

E = N(a + v) + - N{a + V) + \u 2 N{N - l)(a + u)[\n Q - (N - 2)]. (17) 

dE n . 

The extremum condition — = gives a parameter 

d(a + v) 

(a + v) = [l + v 2 (N - l)[ln (?) — (N — 2)]]" 1 / 2 (18) 



at which minimum energy of E is 



E v = N[l + v\N - l)[ln Q - (JV - 2)]] 1 / 2 . (19) 

As established Wu M near the bosonic limit of ground state energy of anyons in 2D parabolic 
well on the factor u, i.e. v — > 0, ground state energy when centre mass motion not fixed has 
a following expression ( here and below in this section we return into normal units of energy 
and length): 

E « %uj [N + N(N — l)v/2] (20) 

for the arbitrary number of particles N in system. So, at v — > we expand expression for 
the energy E v (|l^) on the powers v 2 and by equating term proportional to v 2 of expansion 
to second term in square brackets of (GO), we find that 



^expP+^-^k (21) 



V 



and finally analytical expression for the variational ground state energy is: 

E v = huj N[l + u(N - 1)] 1/2 . (22) 

At the Table 1 we bring exact energies for fermions in 2D parabolic well constructed by 
filling lowest energetic levels of 2D oscillator and energies calculated by formula (p2|). From 
this table one can see that deviation of variational energy E v from exact energy E^xact is 
no more than 6%. 
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TABLES 

TABLE I. The ground state energy of N fermions in 2D parabolic well. 
E ex act - the energies given by exact analytical formula from |14j (in %tu units), 
E v - the energies (in %lo units) calculated by formula (p2[). 
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